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\ field.  The  electric  field  of  the  current  sheet  is  pointing  normally 
toward  the  midplane  of  the  sheet.  This  zeroth-order  perpendicular 
electric  field  is  identified  as  the  primary  electric  field  which  gives 
rise  to  a few  keV  potential  drop  along  auroral  field  lines. 


1.  Introduction 


The  zeroth-order  electric  field  associated  with  a field-aligned 
electron  current  sheet  is  of  fundamental  importance  in  the  formation  of 
discrete  auroras  as  has  been  proposed  by  Kan  and  Akasofu  (1979).  This 
electric  field  is  perpendicular  and  pointing  toward  the  midplane  of  the 
sheet.  Intense  electric  fields  of  this  type  have  been  observed  on 
auroral  field  lines  at  altitudes  of  a few  earth  radii  (R^)  by  Wescott 
et  al.  (1976)  and  Mozer  et  al.  (1977). 

Recently  Kan,  Lee  and  Akasofu  (1979)  showed  that  the  perpendicular 
electric  field  of  the  field-aligned  current  sheet  develops  into  a V-shaped 
potential  structure  due  to  the  high  conductivity  in  the  ionosphere.  As 
a result,  current-carrying  electrons  are  accelerated  downward  along 
magnetic  field  lines  by  the  parallel  electric  field  in  the  V-potential 
region  and  produce  discrete  auroras. 

The  purpose  of  this  paper  is  to  show  the  existence  of  non-neutral 
electron  current  sheets  for  adiabatic  as  well  as  nonadiabatic  electrons 
in  a Maxwellian  plasma  and  to  relate  the  characteristics  of  the  current 
sheet  to  the  thin  auroral  arcs  of  a few  hundred  meters  in  latitudinal 
thickness  (Maggs  and  Davis,  1968).  Non-neutral  electron  current  sheet 
equilibrium  in  a monoenergetic  plasma  has  been  demonstrated  under  the 
assumption  that  the  current  is  exactly  parallel  to  the  external  field 
lines  (Kan  and  Akasofu,  1979).  It  will  be  shown  in  this  paper  that  the 
confinement  of  the  non-neutral  electron  current  sheet  in  a finite  temperature 
Maxwellian  plasma  depends  on  the  diamagnetic  current  perpendicular  to 
the  external  magnetic  field  rather  than  the  parallel  current. 


3 


1 


2.  Basic  Equations 

Equilibrium  configurations  of  exact  charge-neutral  current  sheets 

in  a Maxwellian  plasma  have  been  studied  by  Harris  (1962),  Nicholson 

t 

(1963),  and  Kan  (1972).  However,  a current  sheet  need  not  be  exactly 
neutral,  especially  when  the  current  is  carried  primarily  by  the  electrons. 

A one-dimensional  model  of  the  non-neutral  current  sheet  is  formulated 
on  the  basis  of  the  Vlasov-Maxwell  equations.  The  current  sheet  is 
assumed  parallel  to  the  y-z  plane  and  centered  at  x = 0.  The  external 
magnetic  field  lies  in  the  positive  z-direction.  All  physical  quantities 
may  depend  on  x only. 

The  Vlasov  equation  is  satisfied  by  any  function  of  the  constants 
of  motion  (Longmire,  1963).  However,  there  are  constraints  on  the 
function  and  its  argument  imposed  by  the  prescribed  boundary  conditions 
of  a given  current  sheet  configuration.  In  this  one-dimensional  geometry 
the  constants  of  motion  for  the  s-th  species  are 

H$  “ 2 ™s  <vx  + vy  + vz>  + qs  *(x) 

Pys  = Vy  + (qs/c)  Ay(x)  CD 

Pzs  = Vz  + <qs/c>  Vx) 

The  distribution  function  is  assumed  Maxwellian  and  can  be  written  as 
[ fs(x,#)  = K$  exp  (-cs/6s)  (2) 
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where  the  argument  of  the  exponential  function,  es,  is  an  energy  in 
terms  of  the  constants  of  the  motion,  0$  is  the  thermal  energy,  and 
K$  is  a normalization  constant.  In  general,  can  be  any  linear  combination 
of  H$,  Pys,  Pyz,  P‘s  and  P‘s.  However,  a particular  combination  will  be 
chosen  such  that  the  solutions  satisfy  the  desired  boundary  conditions. 


Nonadiabatic-electron  current  sheet: 

The  field-aligned  current  velocity  V is  assumed  to  be  constant  and  in 
the  positive  z direction.  A diamagnetic  current  is  allowed  to  flow  in 
the  y direction.  Under  these  assumptions,  the  e quantity  for  the  current- 
carrying  electrons  is  chosen  to  be 


ec  = Hc  - vpzc  + 7 mv2  + bPyC/2m 


where  b is  a dimensionless  constant  pertaining  to  the  diamagnetic  current  and 
m is  the  electron  mass.  Upon  substituting  (3)  into  (2)  the  distribution 
function  of  the  current-carrying  electrons  can  be  written  as 


+ bl1* 


.1  + b)  ‘ 
oc 


V,J 


\ ex4vx  * (:  +6)  (vy  * rfb  it  V 


where 


exp[[^  Az  + e U - ♦ ) - pf 


*y^/ec] 


and  ^ $ (x  = 0). 

The  background  electrons  and  background  ions  are  assumed  stationary. 
Density  gradients  are  created  in  the  background  electrons  and  ions  by 
the  electric  field  of  the  current  sheet.  These  density  gradients  produce 
no  current.  The  distribution  function  of  the  background  electrons  may 
be  written  as 

'•  ' "*  'XP  C'  (vx  * vl  * vl)/voel  <«> 

oe 

where 

ne  = Ne  exp  (e*/0e)  (7) 

The  background  ion  distribution  function  can  be  written  similar  to  (6)  as 

f*  ' ;ArT  "1  'XP  L-  <VX  + * vf)/vof  J (8) 

oi 

where 

ni  = Ni  exp  (-e<b/6  - ) (9) 

It  is  convenient  to  introduce  the  following  normalized  quantities 


* w * v„.  * y0 

x«  = x/Ajj,  V*  = V/v0,  E*  = E/E0  (10) 

«/B0,  4*  = 4/B0X„,  ♦*  = ♦/♦„ 
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For  the  auroral  field-aligned  current  sheet  0 * 10"4  and  y * 10-2 , 


the  right  hand  side  of  (12)  and  (13)  can  be  set  to  zero.  One  obtains 


AJ  = B*x*,  A*  = 0 


(17) 


The  number  density  in  (14)  can  be  rewritten  as 


n*  = NJ  exp  ($*  - ♦JJ/eJ  * x*2/2*2 


(18) 


* 2 *2  1/2 

where  SL*  ~ (pe^ya  B ) is  the  normalized  half-width  of  the  nonadiabatic- 


electron  current  sheet. 

Equation  (11)  can  be  solved  numerically  subject  to  the  boundary 
conditions 


♦*  (x*  = o)  = *;  = -** 


(x*  - 0)  ■ 0 


dx 


VS,-  * 0 


(19) 


Lim 


8r  («*)  * o 


where  0*  is  a constant  potential  to  be  determined.  There  are  six 
parameters  in  our  problem  N*  = N*  = N*;  N*;  0*;  9$;  0*;  and  £*.  For  a 
given  set  of  parameters  we  must  find  the  value  of  <frj.  This  is  accomplished 
by  integrating  equation  (11)  starting  from  x*  = 0 using  a trial  value 
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for  <5*.  The  trial  value  for  is  varied  until  the  boundary  conditions 
at  large  x*  in  (19)  are  satisfied. 

Adiabatic-electron  current  sheet: 

Adiabatic  electrons  follow  the  £ x J§  drift  motion  in  the  y direction. 

This  motion  cannot  affect  the  number  density  which  can  only  depend  on  x.  Hence, 
the  background  electrons  can  be  assumed  uniform,  i.e. 

n e = (20) 

Similarly,  the  number  density  of  the  current-carrying  electrons  is 
independent  of  the  electric  field  and  is  determined  externally  by  the  current 
source.  Hence,  it  can  be  assumed  to  have  a Gaussian-shaped  profile,  i.e., 

n*  = NJ  exp  (-x*2/2*2)  (21) 

where  SL*  is  the  half-width  of  the  density  profile.  In  this  case,  the 
thickness  of  the  current  sheet  21*  is  determined  by  the  current  source  as 
long  as  it  is  much  greater  than  the  electron  gyrodiameter.  The  ion  number 
density  is  still  given  by  (16)  as  long  as  the  ions  are  not  completely  magnetized. 

3.  Numerical  Results 

In  this  section  we  show  the  numerical  results  pertaining  to  the  two 
types  of  current  sheets  formulated  in  the  previous  section.  In  Section 
4 we  apply  these  results  to  discrete  auroral  arcs  characterized  by  a 
specific  choice  of  NQ  and  eQ. 
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For  all  examples  in  this  section  we  fix  0*  = 0*  = 0.1;  0*  = 1.0; 

and  N*  = N*  = 1.0  for  the  plasma  parameters  along  auroral  field  lines.  The 

value  of  is  then  determined  by  our  choice  of  2*  and  N*  We  consider 

current  sheets  where  the  half-width  is  on  the  order  of  the  Debye  length, 

i.e. , 2 -v  Aq.  Under  the  quasi-neutrality  approximation,  equation  (11) 

evaluated  at  * = 0 can  be  written  as  <p*  ^ 0?  2n  N*/N*  which  shows  that 

m l ci’ 

the  electric  potential  is  on  the  order  of  the  ion  thermal  energy.  For  a 

given  2*  we  choose  M*  such  that  <}>*  is  on  the  order  of  0?. 

c m l 

Nonadiabatic-electron  current  sheet: 

To  illustrate  nonadiabatic-electron  current  sheet  solutions  two 

examples  are  considered.  The  first  example  is  present  in  Figure  1 with 

2*  = 0.5  and  N*  = 4.06.  Figure  la  shows  the  number  densities  of  the 

ions  and  electrons.  The  width  of  the  electron  density  profile  is  on  the 

order  of  the  Oebye  length.  Note  that  the  current-carrying  electrons  are 

confined  within  x*  = 1.5  while  the  background  electrons  are  pushed  out 

of  the  current  sheet  by  the  electric  potential  $*.  This  behavior  is  the 

cause  of  the  dip  in  the  electron  density  profile  near  x*  = 1.8.  The 

ions  are  attracted  to  the  electron  current  sheet  and  partially  shield 

the  electric  field.  As  shown  in  Figure  lb,  the  depth  of  the  potential 

profile  is  on  the  order  of  the  ion  thermal  energy  (ra  = 1.2  ©•)•  The 

maximum  electric  field  is  E?  = 0.8. 

m 

The  second  example  is  shown  in  Figure  2 with  2*  = 0.2  and 
N*  = 7.7.  The  depth  of  the  potential  well  <|>*  is  about  1.2.  In  this 
case,  the  current-carrying  electrons  have  a higher  number  density  and  the 
current  sheet  is  confined  to  a smaller  region  when  compared  to  the  example 
shown  in  Figure  1.  As  a result,  the  current  sheet  supports  a larger 
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electric  field,  E*  = 1.4.  Note  that  the  dip  in  the  electron  density  profile 
in  Figure  2a  is  deeper  than  the  dip  in  Figure  la. 


Adiabatic-electron  current  sheet: 

Figure  3 shows  an  adiabatic-electron  current  sheet  solution  with  l* 

= 1.0  and  N*  = 3.3.  The  choice  of  N*  = 3.3  gives  <)>*  = 1.2.  The  width 
of  the  electron  density  profile  is  on  the  order  of  the  Debye  length. 

The  maximum  electric  field  is  E*  = 0.6.  There  is  no  dip  in  the  electron 
density  profile  because  the  electron  distribution  functions  are  not 
affected  by  the  formation  of  the  electric  potential. 

4.  Application  to  Discrete  Auroral  Arcs 

The  numerical  solutions  we  have  obtained  are  valid  for  any  choice 

of  N.  and  0„.  ti 

o o 

Consider  an  auroral  field-aligned  electron  current  sheet  at  1 Rg 
altitude  where  NQ  = 1 cm  and  0Q  = 5 kev  are  the  particle  number 
density  and  the  ion  thermal  energy,  respectively.  The  magnetic  field 
at  1 is  about  0.06  Gauss.  The  Debye  length  is  Aq  ^ 500  m and  the  electron 
gyroradius  is  pe  ^ 12  m in  the  absence  of  the  electric  field.  In  this 
case,  electron  motions  are  adiabatic.  Hence,  we  can  apply  the 
solution  shown  in  Figure  3 to  obtain  0C  = 0g  = 0.5  kev,  0.  = 5 kev, 

-3  .3 

Ng  = N-  = 1 cm  , and  Nc  = 3.3  cm  . The  half-width  of  the  current  layer 
is  * 500  m.  The  depth  of  the  potential  profile  is  ^ 6 kev  and  the  maximum 
electric  field  is  ^ 6 V/m. 

As  a second  example,  we  consider  a current  sheet  at  2 R^  altitude 

_3 

in  which  NQ  = 0.3  cm  , 0Q  = 5 kev  are  the  typical  particle  density 
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and  the  typical  ion  thermal  energy,  respectively.  The  magnetic  field  at 
2 Rj:  is  ~ 0.02  Gauss.  In  this  case,  \q  * 860  m and  pg  ~ 45  m in  the 
absence  of  the  electric  field.  It  can  be  shown  from  numerical  calculation 
of  particle  trajectories  that  the  electrons  are  nonadiabatic  if  the  half- 
width £ is  less  than  ~ 200  m.  The  nonadiabatic-electron  current  sheet 
solution  shown  in  Figure  2 is  applicable  since  £ = 0.2  Ag  * 170  m.  For 
the  solution  in  Figure  2 we  have  Ng  = N-  = 0.3  cnf^;  Nc  = 2.3  cm"'*; 

9e  = ®c  s kev»  and  9i  = ^ lcev*  ^he  dePtfl  °f  tfie  potential  profile 
is  ^ 6 kev  and  the  maximum  electric  field  in  the  current  sheet  is  * 8 V/m. 
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5.  Conclusion 

We  have  shown  the  existence  of  non-neutral  field-aligned  electron 
current  sheets  in  a Maxwellian  plasma  for  adiabatic  as  well  as  nonadiabatic 
electrons.  The  electric  fields  are  confined  within  the  current 
sheet  due  to  the  shielding  effect  of  the  background  ions.  If  the  current- 
carrying  electrons  are  nonadiabatic,  the  electron  number  density  profile 
depends  on  the  electrostatic  potential.  On  the  other  hand,  if  the  electrons 
are  adiabatic,  the  electron  number  density  profile  is  determined  externally 
by  the  current  source. 

The  results  show  that  the  total  potential  energy  in  the  non-neutral 
current  sheet  is  on  the  order  of  the  ion  thermal  energy  which  ranges  from 
1 to  10  keV  in  the  plasma  sheet.  The  model  gives  a maximum  electric  field 
ranging  from  1 to  10  V/m.  Due  to  the  short-circuiting  effect  of  the 
conductive  ionosphere,  the  total  potential  of  the  current  sheet  will  be 
forced  to  be  distributed  along  magnetic  field  lines  as  has  been  shown  by 
Kan  et  al.  (1979).  This  energy  is  comparable  to  the  energy  of  the 
precipitating  auroral  electrons. 
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Figure  Captions 


Figure  1 


Figure  2 


A nonadiabatic-electron  current  sheet  solution  for 
2/A.q  = 0.5. 

(a)  The  dimensionless  electron  and  ion  number  density 
profiles  as  a function  of  the  dimensionless  distance, 
x/Xg,  from  the  center  of  the  current  sheet. 

(b)  The  electric  potential  and  perpendicular  electric 

field  profiles  as  a function  of  x/Xg.  For  the  conditions 

.0 

at  “n.  2 Rg  altitude  on  auroral  field  lines:  NQ  = 0.3  cm 
and  0O  = 5 keV  giving  Xg  = 860  m;  0O  = 5 kV;  and  EQ  = 5.5  V/m. 

A nonadiabatic-electron  current  sheet  for  £/Xg  = 0.2. 

(a)  The  dimensionless  ion  and  electron  number  density  profiles 

as  a function  of  x/Xg.  The  number  density  of  current  carrying 
electrons  is  greater  than  for  Figure  1,  and  the  current 
carrying  electrons  are  confined  to  a smaller  region. 

(b)  The  electric  potential  and  perpendicular  electric  field 
as  a function  of  x/Xg.  The  electric  field  supported  by 
the  current  sheet  is  larger  than  in  Figure  1.  For  the 
condition  at  * 2 Rg  altitude  on  auroral  field  lines: 

NQ  = 0.3  cm”'*  and  0Q  = 5 keV  giving  Xg  = 860  m;  <j>Q  = 5 kv; 
and  EQ  s 5.5  V/m. 
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Figure  3 An  adiabatic-electron  current  sheet  solution  for  2/A.g  = 1.0. 

(a)  The  dimensionless  ion  and  electron  number  density  profiles 
as  a function  of  x/\q. 

(b)  The  dimensionless  electric  potential  and  perpendicular 

electric  field  as  a function  of  x/\g.  For  the  conditions 

at  ^ 1 Rc  altitude  on  auroral  field  lines:  N„  = 1 cm’^ 
t o 

and  0O  = 5 keV  giving  Xg  = 500  m;  = 5 kV;  and  EQ  = 9.5  V/m. 
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